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Abstract 

A series of semi-inclusive deep-inelastic scattering measurements on deuterium, he- 
lium, neon, krypton, and xenon targets has been performed in order to study ha- 
dronization. The data were collected with the HERMES detector at the DESY 
laboratory using a 27.6 GeV positron or electron beam. Hadron multiplicities on 
nucleus A relative to those on the deuteron, R\, are presented for various hadrons 
(ir + , 7T _ , 7T°, K + , K~ , p, and p) as a function of the virtual-photon energy v, the 
fraction z of this energy transferred to the hadron, the photon virtuality Q 2 , and 
the hadron transverse momentum squared p 2 . The data reveal a systematic decrease 
of R h A with the mass number A for each hadron type h. Furthermore, R\ increases 
(decreases) with increasing values of v (z), increases slightly with increasing Q 2 , and 
is almost independent of p 2 , except at large values oip 2 . For pions two-dimensional 
distributions also are presented. These indicate that the dependences of R\ on v 
and z can largely be described as a dependence on a single variable L c , which is a 
combination of v and z. The dependence on L c suggests in which kinematic condi- 
tions partonic and hadronic mechanisms may be dominant. The behaviour of R\ at 
large p 2 constitutes tentative evidence for a partonic energy-loss mechanism. The 
A-dependence of R\ is investigated as a function of v, z, and of L c . It approximately 
follows an A a form with a ~ 0.5 — 0.6. 

Key words: nuclei, quarks, hadron production, hadronization, attenuation, 
A-dependence 



1 Introduction 



After decades of extensive study, understanding the confinement of quarks and 
gluons in hadrons still is one of the great challenges in hadronic physics. To 

1 Present address: 36 Mizzen Circle, Hampton, Virginia 23664, USA 
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uncover its nature, hadronic reactions in a nuclear medium, either cold or hot, 
are studied. Typical examples are the measurements of hadron production on 
nuclear targets in semi-inclusive deep-inelastic scattering of leptons [1,2,3,4,5] 
and the jet-quenching and parton energy-loss phenomena observed in ultra- 
relativistic heavy-ion collisions [6,7]. In each case hadron yields are observed 
that are different from those observed in the corresponding reactions on free 
nucleons. 

The process that leads from the partons produced in the elementary inter- 
action to the hadrons observed experimentally is commonly referred to as 
hadronization or fragmentation. According to theoretical estimates the hadro- 
nization process occurs over length scales varying from less than a femtometer 
to several tens of femtometers. At these length scales the magnitude of the 
strong coupling constant is such that perturbative techniques cannot be ap- 
plied. Hence, hadronization is an intrinsically non-perturbative QCD process, 
for which only approximate theoretical approaches are presently available. Ex- 
perimental data are vital for supporting these theoretical developments, since 
they can be used to gauge or guide the calculations. 

The hadronization process in free space has been studied extensively in e + e~ 
annihilation experiments [8]. As a result the spectra of particles produced and 
their kinematic dependences are rather well known. However, little is known 
about the space-time evolution of the process. Semi-inclusive production of 
hadrons in deep-inelastic scattering of leptons from atomic nuclei provides a 
way to investigate this space-time development. Leptoproduction of hadrons 
has the virtue that the energy and the momentum of the struck parton are 
well determined, as they are tagged by the scattered lepton. By using nuclei 
of increasing size one can investigate the time development of hadronization. 
If hadronization occurs quickly, i.e., if the hadrons are produced at small 
distances compared to the size of atomic nuclei, the relevant interactions in 
the nuclear environment involve well-known hadronic cross sections such as 
the ones for pion-nucleon interactions. If, in contrast, hadronization occurs 
over large distances, the relevant interactions are partonic and involve the 
emission of gluons and quark-antiquark pairs. The two mechanisms lead to 
different predictions for the decrease in hadron yield, known as attenuation, 
on nuclei as compared to that on free nucleons. 

Most likely, a combination of these two mechanisms contributes to the ob- 
served attenuation of hadron yields on nuclei. This expectation has led to a 
range of phenomenological approaches, which are briefly summarized in the 
next section. The available calculations cover a range of possible mechanisms 
(partonic energy loss or hadronic absorption) and time (length) scales (from 
less than 1 fm to more than 10 fm), which results in different dependences 
on the various kinematic variables. In order to distinguish between these cal- 
culations precise hadron attenuation data are needed as a function of several 
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kinematic variables for a range of nuclei and for several hadron types. 

Exploratory measurements were first performed at SLAC [1] and later by the 
EMC [2] and E665 [3] collaborations. More recently many more data have 
been collected by the HERMES collaboration at DESY and the CLAS col- 
laboration at the Thomas Jefferson National Accelerator Facility [9]. Some 
of the HERMES data have been published already [4,5]. The CLAS data are 
presently being analysed [10]. In this paper we present the full results from 
HERMES on the multiplicities for the production of pions, kaons, protons, 
and antiprotons on helium, neon, krypton, and xenon targets relative to those 
for deuterium. It goes beyond the scope of the present paper to compare the 
data to the available theoretical calculations. Instead, prominent features of 
the data are identified and used to address two key issues in the study of 
hadronization: (i) what are the time or length scales of the process, and (ii) 
what are the mechanisms that lead to nuclear attenuation? 

The paper is organized as follows. In section 2 the theoretical framework is 
described, and some representative theoretical models are summarized. In sec- 
tion 3 those aspects of the HERMES experiment that are relevant to the 
present measurements are presented. In section 4 the data analysis is dis- 
cussed, including the corrections to the raw data and the systematic uncer- 
tainties. The results for the attenuation as a function of various kinematic 
variables are presented and discussed in section 5. This section has several 
subsections in which the features of the data, especially those related to the 
relevant time scales and mechanisms, and the A-dependence, are discussed 
separately. The results are summarized in the last section, which also lists the 
conclusions. 



2 Theoretical framework 

In order to put the experimental results into perspective, in this section the 
concepts that are used in the study of hadronization are presented, and mod- 
els that have been developed to describe the experimental results are briefly 
discussed. It is emphasized that the latter is meant only to illustrate potential 
interpretations of the data. The experimental data and the features that they 
exhibit are the genuine subject. 

The hadronization process in a nuclear medium can be studied by means of 
semi-inclusive deep-inelastic scattering (SIDIS) of electrons or positrons from 
nuclei. For that purpose the multiplicity ratio R\ is introduced, which is de- 
fined as the ratio of the number of hadrons h produced per deep-inelastic 
scattering (DIS) event on a nuclear target with mass number A to that for a 
deuterium (D) target. Figure 1 illustrates the definition of the relevant lep- 
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Fig. 1. Kinematic planes for hadron production in semi- inclusive deep-inelastic scat- 
tering and definitions of the relevant lepton and hadron variables. The quantities 
k = (E, k) and k' = (E' , k') are the four-momenta of the incident and scattered lep- 
ton, (Eh,Ph) is the four-momentum of the produced hadron, and pt is the transverse 
momentum of the hadron. 

ton and hadron kinematic variables. The ratio R h A depends on the leptonic 
variables u, the energy of the virtual photon, and Q 2 , the negative of the 
four-momentum of the virtual photon squared, and on the hadronic variables 
z = Eh/v, the fraction of the virtual-photon energy carried by the hadron, and 
p 2 , the square of the hadron momentum component transverse to the direction 
of the virtual photon. Thus R\ can be written as: 



R h A (u,Q 2 ,z,p 2 ] 



( N h (u,Q 2 ,z,p 2 ) 
\ N e (v,Q 2 ) 
f N h (u,Q 2 ,z,p 2 ) \ ' 
V N*(u,Q 2 ) J D 



with N h (v, Q 2 , z, p\ ) the number of semi-inclusive hadrons at given (V, Q 2 , z, p\) . 
and N e {v, Q 2 ) the number of inclusive DIS leptons at (u, Q 2 ). Implicit in this 
definition is the integration over the angle between the lepton scattering 
plane and the hadron production plane (see Fig. 1). 

Experiments at large values of v [2,3] give values R A ~ 1.0 within the exper- 
imental uncertainty. This is interpreted as an indication that nuclear effects 
are negligible in that region. At lower values of v the value of R A has been 
found to be well below unity [1,4,5]. 

Even though hadronization is not yet quantitatively understood, it is generally 
assumed that the following processes play a role in leptoproduction of hadrons. 
After a quark in a nucleon has absorbed the virtual photon, it can lose energy 
by scattering from other quarks and by radiating gluons. As a result the value 
of R\ may be influenced. A change in R h A can also result from the quark 
in a nucleus having a different distribution function as a result of partonic 
rescaling. These two effects will be called partonic mechanisms. 

After a certain time a colourless object, called a prehadron, is formed, which 
has the quantum numbers of the final hadron, but not yet its full wave func- 
tion. The concept of colour transparency [11] is closely related to this. The 
prehadron then evolves over some time into the physical hadron. In the Lund 
model [12] the process of (pre)hadron formation is described as the building of 
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a colour field (string) between the struck quark and the residual system. The 
string gets stretched and breaks into smaller pieces, each with an (anti-)quark 
at both ends. The Lund model gives predictions (see, e.g., Ref. [13]) for the 
time t c , and the corresponding (constituent or formation) length l c , that it 
takes for the prehadron to be formed, and the time th required for the final 
physical hadron to be formed!" 2 ! 

The average value of the formation length l c , denoted by L c , is given in the 
Lund model by 

L c = /(*)-, (2) 

K 

where f(z) is a function that goes as 1 — z for large values of z (this behaviour 
is also found in other models, and is due to the fact that in the limit z — > 1 the 
struck quark cannot have lost any energy when producing the final hadron), 
and has a broad maximum of about 0.4 around z = 0.35. The quantity k 
is the string tension, which reflects the energy loss of the leading quark per 
unit length, usually taken as K — 1.0 GeV/fm. Thus, at values of v of 5 — 25 
GeV, the prehadron formation takes place over distances of about 1 — 10 fm, 
comparable to the size of a nucleus. The values of 1^ are larger than l c , so the 
final hadron is often formed outside the nucleus. 

If the (pre)hadron is formed inside the nucleus, it can experience hadronic 
interactions, generally called final-state interactions (FSI). For discussing the 
effects of these, we will not discriminate between prehadrons and hadrons, and 
just talk about hadrons and hadronic effects, although the relevant cross sec- 
tions may be different for hadrons and prehadrons. The first effect of FSI may 
be the rescattering of the hadron from the nucleons in the nucleus, the hadron 
losing energy and possibly 'generating' other, mainly low energy, hadrons (if 
these are mesons, the word generate can be taken literally; if these are nu- 
cleons, generated means that they are now in a continuum state). The effect 
is thus a loss of hadrons at a given value of z, and an increase of the same 
hadron, but possibly also of other hadrons, at lower z. Another effect of FSI 
can be that the original hadron is absorbed. This is usually accompanied by 
the emission of other, again mainly low-energy, hadrons. So also in this case 
there is a loss of hadrons at given z, and an increase of other hadrons at lower 
z. For a full description of all these effects, coupled- channel calculations should 
be performed. However, since these are rather complicated, in models usually 
only absorption is considered. Clearly this neglects the generation of hadrons 
at lower values of z. 



2 In the literature the name 'formation time or length' has sometimes different 
meanings and also different symbols are used. Here the name 'formation length' (l c ) 
is used for what is often called the prehadron (or constituent) formation length, and 
'hadron formation length' (Z/J for the hadron (or yo-yo) formation length (see, e.g., 
Refs. [13,14]). 
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The change from partonic mechanisms to hadronic mechanisms, which de- 
pends on the formation length relative to the size of the nucleus, influences 
R h A . Thus the use of various nuclear targets allows one to investigate the space- 
time development of the hadronization process. 

At present, reliable QCD calculations of quark hadronization (fragmenta- 
tion) can not yet be performed because of the major role of 'soft', i.e., non- 
perturbative processes. For that reason various types of models have been 
developed. 

Phenomenological models [15,16,17,18] use formation times/lengths and ab- 
sorption cross sections for the various hadrons in the nuclear medium. Various 
formulae for the formation lengths have been used, and in the more advanced 
versions two length scales, l c and lh, are distinguished, as well as different ab- 
sorption cross sections for the quark, the prehadron, and the final hadron. The 
absorption cross sections are usually adjusted to obtain the best description 
of the v- and z-dependence of the experimental data. 

Other (QCD-inspired) models focus on the energy loss that the struck quark 
experiences in the nuclear environment [19,20,21,22]. In Refs. [20,21] twist-4 
contributions to the fragmentation functions resulting from multiple scatte- 
ring and gluon bremsstrahlung in a nuclear medium are calculated. A nuclear 
attenuation proportional to A 2 / 3 is predicted, where the power 2/3 results 
from the coherence of the gluon radiation process [23], which gives an induced 
radiative energy loss of a quark traversing a length L of matter proportional 
to L 2 [24] . No hadron absorption is included, as it is assumed that the hadron 
is formed outside the nucleus. 

The effect of a finite formation length is included in Ref. [22], in which fragmen- 
tation functions are calculated that are again modified due to partonic energy 
loss in the nuclear medium during the time t c . By using the quark energy loss 
determined from Drell-Yan data, a reasonable agreement with existing data is 
found. In order to keep the approach as simple as possible, absorption of the 
produced (pre)hadron is not taken into account. 

Another class of models includes (pre) hadron absorption, with or without a 
description of what happens in the hadronization process. In Ref. [25] R\ is 
described in terms of medium-modified fragmentation functions supplemented 
by nuclear absorption. A parton-rescaling model that has also been used to 
describe the EMC effect is used to describe the nuclear modification of the 
fragmentation functions. The (average) formation length is taken from the 
Lund model [12]. Ref. [26] calculates the nuclear attenuation of the (leading) 
hadron with z > 0.5 by including as major ingredients the formation length l c 
and an absorption cross section of the prehadron. The effect of quark energy 
loss is found to be small. In Ref. [27] the hadron attenuation is investigated 
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within the framework of the Boltzmann-Uehling-Uhlenbeck transport model. 
Since this is a coupled-channel approach, hadrons are not only absorbed, but 
can also be produced. Some choices for the formation time t c , including taking 
it to be zero, are studied. 

The theoretical calculations have been compared to the data from Refs. [2,4,5]. 
Notwithstanding their different and sometimes orthogonal ingredients, all mod- 
els reproduce the global features of the data. In the case of the ^-dependence, 
which is best described, the reason may be that the decreasing attenuation 
with increasing value of v is largely due to a simple increase of the formation 
time t c with v in the target rest frame due to Lorentz dilatation. The depen- 
dence of R\ on z and A, which in general is less well described, may be more 
discriminating, especially when more detailed data are available. 



3 Experiment 

The measurements were performed with the HERMES spectrometer [28] us- 
ing a 27.6 GeV positron or electron beam stored in the HERA ring at DESY. 
Some data were collected using a 12.0 GeV positron beam [29,30], but since 
the amount of data was much less and only pions or all hadrons were iden- 
tified, they are not included here. Typical beam currents were 40 mA down 
to 5 mA. The spectrometer consists of two mirror-symmetric halves, located 
above and below the lepton beam pipe. A flux-exclusion plate in the midplane 
of the magnet prevents deflection of the lepton (and proton) beams passing 
through the center of the spectrometer. The scattered leptons and the pro- 
duced hadrons were detected within an angular acceptance of ± 170 mrad 
horizontally and ± (40 - 140) mrad vertically. The lepton trigger was formed 
by a coincidence between signals from three scintillator hodoscope planes and 
a lead-glass calorimeter. A minimum energy deposit in the latter of 3.5 GeV 
(1.4 GeV) for unpolarized (polarized) target runs was required. 

The data were collected during the years 1999, 2000, 2004, and 2005, using 
unpolarized nuclear (He, Ne, Kr, Xe) and polarized or unpolarized deuterium 
(D) gaseous targets internal to the storage ring (see Table 1). In 1997 also data 
on nitrogen were taken, but since at that time the RICH detector (see below) 
was not yet available, only data for pions and all hadrons together could be 
presented [4]. Therefore these data are not included here. The yields from 
polarized deuterium were summed over the two spin orientations. The target 
gases were injected into a 40 cm long tubular open-ended storage cell. Using 
an Unpolarized Gas Feed System [31] it is possible to provide D, He, N, Ne, 
Kr, or Xe targets with relatively high areal densities (between 10 14 and 10 17 
nucleons/cm 2 ), resulting in luminosities ranging from 10 31 to 10 33 cm~ 2 s _1 . 
Such high-density runs were taken at the end of HERA fills, with typical 
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currents of 15 to 5 rriA and beam lifetimes of one hour. This made it possible 
to accumulate the data for these targets in only a few days of integrated 
beam time. The luminosity was measured using elastic scattering of the beam 
particles from the electrons in the target gas: Bhabha scattering for a positron 
beam, M0ller scattering for an electron beam [32]. There are several particle 

Table 1 

Overview of the HERMES nuclear attenuation measurements. 



Year 


E (GeV) 


Target 


Identified hadrons 


Ref. 


1997 


27.6 


D, N 


h ± ,7T ± 


[4], [5] 


1999 


27.6 


D, Kr 


/l ± , 7T , 7T°, K^, p, p 


[5], this work 


2000 


27.6 


D, He, Ne 


ir^jtr , K^, p, p 


this work 


2000 


12.0 


D, N, Kr 


h ± , ^ 


[29] 


2004 


27.6 


D, Kr, Xe 


ir^, if ^ 7r°, p, p 


this work 


2005 


27.6 


D, Kr, Xe 


ir^, K^jTr , p, p 


this work 



identification (PID) detectors in the HERMES spectrometer. Details on the 
performance and use of these PID detectors can be found in Ref. [33] . Electrons 
and positrons are identified by combining the information from a lead-glass 
calorimeter, a scintillator hodoscope preceded by two radiation lengths of lead 
(the pre-shower detector), and a transition-radiation detector. 

The identification of charged pions, kaons, protons, and antiprotons is ac- 
complished using the information from the Ring-Imaging Cerenkov detector 
(RICH) [34], which replaces a threshold Cerenkov counter used in the pre- 
viously reported measurements on nitrogen [4]. The RICH detector uses two 
radiators, a 5 cm thick wall of silica-aerogel tiles followed by a large volume of 
C4F10 gas, to provide separation of pions, kaons, and (anti)protons. Together 
these provide good particle identification for charged hadrons in the momen- 
tum range 2 < p < 15 GeV, with limited contamination from misidentified 
hadrons. The identification efficiencies and contaminations for pions, kaons, 
protons, and antiprotons were determined in a Monte Carlo simulation as 
a function of the hadron momentum and multiplicity in the relevant detec- 
tor half. These performance parameters were verified in a limited kinematic 
domain using known particle species from identified resonance decays. They 
were used in a matrix method to unfold the true hadron distributions from 
the measured ones. Systematic uncertainties in the unfolding were estimated 
by using matrices determined in different ways, see Refs. [34,35] for details. 

The electromagnetic calorimeter [36] provides neutral pion identification through 
the detection of two clusters without an accompanying track, originating from 
the two decay photons. 
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4 Data Analysis 



The analysis procedure is similar to the one described in detail in Refs. [4,5], 
where the nitrogen and first krypton data were presented. Since the publication 
of Ref. [5], more data on krypton were taken and all data were analysed in a 
wider kinematic range. 

The hadron multiplicity ratio R\ as defined in Eq. 1 was determined as a 
function of the leptonic (Q 2 and v) and hadronic (z and p 2 ) variables for all 
identified particles and all targets. The kinematic constraints imposed on the 
scattered leptons were identical for all analysed data: Q 2 > 1 GeV 2 , W = 
y/2Mu + M' 2 -Q' 2 > 2 GeV (where M is the nucleon mass) for the invariant 
mass of the photon-nucleon system, and y = vjE < 0.85 for the energy 
fraction of the virtual photon. The constraints on W and y were applied in 
order to exclude nucleon resonances and to limit the magnitude of the radiative 
corrections to R\, respectively. The resulting value of xbj = Q 2 jlMv ranged 
from 0.023 to 0.8. 

As mentioned in the previous section, charged hadrons were identified in the 
momentum range 2.0 — 15.0 GeV by using the RICH detector. For the iden- 
tification of the neutral pions through the detection of their decay photons, 
each of the two photon clusters was required to have an energy of at least 
0.8 GeV. The background was evaluated in each kinematic bin by fitting the 
two-photon invariant mass spectrum with a Gaussian plus a polynomial that 
fits the shape of the background due to uncorrelated photons. The number 
of detected neutral pions was obtained by integrating the peak, corrected for 
background, over a ±2<r range with respect to the centroid of the Gaussian. 
The low momentum limit for the neutral pions was set at 2.5 GeV in order to 
reduce backgrounds. 

The integrated luminosities for all years and targets are listed in Table 2. 
Typical numbers of observed DIS leptons and identified hadrons are listed in 
Table 3. 

Most of the systematic uncertainties related to the detector, the reconstruc- 
tion efficiencies and particle identification practically cancel in the ratio of the 
multiplicities. In determining the multiplicity ratios, deuterium data collected 
in the same year as the data for the heavier target were used to avoid uncer- 
tainties due to possible different conditions or functioning of the HERMES 
spectrometer during the years. It was verified that the multiplicity ratios ob- 
tained in different years were consistent within the statistical and systematic 
uncertainties. 

The multiplicity ratios were also inspected as a function of the hadron angles 
6 X and 9 y with respect to the beam direction in order to investigate whether 
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Table 2 

Integrated luminosities (in pb" 1 ) for the various data sets. 



Target 


1999 


2000 


2004 


2005 


Sum 


D 


32.3 


119.7 


35.7 


61.7 


249.4 


He 




27.9 






27.9 


Ne 




84.2 






84.2 


Kr 


26.1 




29.5 


21.1 


76.7 


Xe 






21.2 


21.4 


42.6 



Table 3 

Number of DIS leptons and identified hadrons collected on deuterium and kryp- 
ton targets in 1999, 2004, and 2005 combined. The numbers are for the following 
kinematic constraints: Q 2 > 1 GeV 2 , v > 6 GeV, W > 2 GeV, and z > 0.2. 



Target 


DIS 


7T+ 


7T 


7T° 


K + 


K~ 


V 


V 


D 


6669k 


706k 


575k 


232k 


146k 


62k 


131k 


23k 


Kr 


3516k 


286k 


232k 


90k 


68k 


26k 


69k 


8k 



the values of R A depend on the geometrical acceptance of the spectrometer. 
After applying (small) corrections for changes in average kinematics with 9 X 
or 9 y , no effect was found beyond the statistical and systematic uncertainties. 
The dependence of the value of R\ on the azimuthal angle <j) of the hadron (see 
Fig. 1) was also investigated, since it is known that the SIDIS cross section on 
the proton and deuteron depends on this variable. It was found that R A was 
constant as a function of within the statistical and systematic uncertainties. 

The data for the multiplicity ratios were corrected for radiative processes in 
the manner described in Ref. [37]. The code of Ref. [38] was modified to 
include the measured SIDIS cross sections. The radiative corrections (RC) 
were applied to both the inclusive and the semi- inclusive parts in Eq. 1. For 
the inclusive cross sections elastic, quasi-elastic, and inelastic processes need 
to be taken into account, whereas for the semi- inclusive ones only inelastic 
radiative processes contribute. The correction for the ratio of the latter was 
taken to be independent of z. Since the inelastic radiative effects are almost 
the same for the nuclei A and D, the size of the radiative corrections applied 
to R\ was small over most of the kinematic range. Only in kinematic regions 
of DIS where the elastic and quasi-elastic tails are non- negligible, i.e., at the 
highest value of v and lowest value of Q 2 (low XBj), is there a noticeable effect 
on R\, with a maximum (increase) of R\ of about 7% for xenon and krypton, 
4.5% for neon, and 1% for helium. 
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Since the usual interpretation of the definition of R\ (see Eq. 1) is that it 
should only include hadrons formed in the fragmentation process, a correction 
has to be made for measured hadrons that are the decay products of directly 
produced mesons. The main effect is on the charged-pion multiplicities as a 
result of the decay of exclusively produced p° mesons (for pions from other 
mesons and for other hadrons the contribution is small). That may affect the 
multiplicities for positive (negative) pions by an amount ranging from about 
1% at low z up to 30% (45%) at high values of z , as estimated from a Monte 
Carlo simulation. The effect on the multiplicity ratio R\ is much smaller, 
but does not cancel completely since the p° mesons also interact with the 
nuclear medium. Taking into account the measured nuclear transparency [39] 
for p° mesons, the maximum remaining effect on R\, which occurs for ^-values 
of 0.7-0.8, was estimated to be about 2(4)%, 3(5)%, 3.5(6)%, and 4(7)% in 
the case of helium, neon, krypton, and xenon, respectively. The first (second) 
number applies for positive (negative) pions. These values were included in 
the systematic uncertainties. 

The total systematic uncertainties include the uncertainties of hadron identi- 
fication (1.5% for neutral pions, 0.5% for charged pions, 2% for kaons, 2% for 
protons, and 6% for antiprotons) , overall efficiency (< 2%), p°-meson produc- 
tion for positive (0.3% - 4%) and negative (0.3% - 7%) pions, and the effects of 
using different parameterizations of fragmentation functions and distribution 
functions in the RC calculations (< 2%). 



5 Results and Discussion 

In this section the experimental results are presented and the dependences of 
the multiplicity ratios R\ on the various kinematic variables and the mass 
number A of the nucleus are discussed. Unless specified otherwise the data 
are shown with the following constraints: v > 6.0 GeV, z > 0.2, and xf > 0, 
where xp is given by 

XF =p|/ p ~, (3) 

with pjj being the component of the hadron momentum in the direction of 
the momentum transfer in the virtual-photon nucleon center-of-mass system. 
Together with that on z, the constraint onip will reduce possible contributions 
from target fragmentation. 

5. 1 Multiplicity ratio for identified hadrons 

Figures 2-5 show the dependence of R\ on u, z, Q 2 , and p 2 for the various 
nuclei for all identified hadrons: positively charged (pions, kaons, and protons), 
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negatively charged (pions, kaons, and anti-protons), and neutral ones (pions). 
The inner error bars in these figures represent the statistical uncertainties, 
while the outer ones are for the total uncertainty (statistical plus systematic, 
added quadratically). The systematic uncertainty is mainly a scale uncertainty, 
affecting the values of R\ for the various values of u, z, Q 2 , or p 2 in the same 
way. 

In presenting the results for R\ as functions of one of the four independent 
variables (v, z, Q 2 , p 2 ) only, R\ was integrated (within the acceptance of 
the experiment) over the others. Because of acceptance effects and because 
in general the dependence of R\ on u, z, Q 2 , and p\ does not factorize 
(R h A (y,z,Q 2 ,p 2 ) ^ R 1 (u)R 2 (z)R 3 (Q 2 )R 4 (p 2 )) J this integration may introduce 
false dependences. This is mainly relevant in case of v and z, where the av- 
erage value of v (z) changes non- negligibly depending on the value of z (u). 
Table A.l in the Appendix gives an indication of the size of these correlations 
by listing the average values of the kinematic quantities that were integrated 
over for the various dependences in the case of produced pions. All data are 
available in detail from Ref. [40]. 

Before discussing in the following subsections the dependence of R h A on the 
kinematic variables z, Q 2 , and pf in detail, first some global features of the 
data are discussed. 

The basic feature of the data is the decrease of R\ with increasing value of 
the mass number A of the nucleus. Qualitatively this is understood as being 
due to increased partonic (quark energy loss) or hadronic (absorption) effects. 
Furthermore, there is a large similarity between the data for ir + , ir~ , and it , 
and a clear difference between those for K + and K~ , and those for p and p. 
Also here there are some simple arguments to explain these features at least 
qualitatively. 

Since we use (almost) isoscalar targets and the production of ir + and n~ 
on protons or neutrons is only slightly different, both the production and 
absorption of pions is in the first instance independent of their charge. The 
values of R\ for K + and K~ show a similar behaviour as a function of the 
various variables, but R A is almost everywhere smaller than R A + . A positive 
kaon can be produced directly from the struck quark (in the language of 
string breaking models it is a rank 1 hadron), but a negative kaon can only be 
produced in more complicated string breakings (rank 2 or higher), except at 
small values of XBj, where sea quarks start to play a larger role. This is reflected 
in the production rate on deuterium being much larger for K + than for K~ , 
see Table 3, and leads to a steeperr dependence of the K~ fragmentation 
function on z, and hence a reduced production if the parton has lost energy 
before hadronization. Also, due to their quark content, nuclear absorption 
cross sections are larger for K~ than for K + . Thus both parton energy loss 
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Fig. 2. Values of R\ for positively charged hadrons as a function of z/, z, and Q 2 . 
The data as a function of ^ are shown for v > 4 GeV and those as a function of 
z for z > 0.1. The inner error bars represent the statistical uncertainty, while the 
outer ones show the total uncertainty. 



and absorption of the produced kaon can qualitatively explain the observed 
difference between + and R\ . However, when comparing the multiplicity 
ratios for pions and kaons, it is seen that R\ wflj, whereas R% > R* A . 
Given that both pions and K + particles are rank-1, and K~ rank 2 or higher, 
and that nuclear absorption cross sections for both K + and K~ are smaller 
than for pions, these features are not readily explained by the behaviour of 
fragmentation functions or absorption cross sections. 

The results for protons cannot really be related to those for any of the other 
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Fig. 3. Values of R\ for negatively charged hadrons as a function of u, z, and Q 2 . 
The data as a function of v are shown for v > 4 GeV and those as a function of z 
for z > 0.1. Error bars as in Fig. 2. 

particles. Because protons are already present in a nucleus, an appreciable 
fraction of them may not come from hadronization. This is reflected in the 
very large difference in production of p and p on deuterium, see Table 3. 
Furthermore, as discussed in section 2, in final-state interactions they generally 
are not absorbed, but give rise to more nucleons (both protons and neutrons), 
thus possibly even increasing R\ at lower z. 

Antiprotons feature a rather strong attenuation, which might be attributed to 
the relatively large pN cross section. 
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Fig. 4. Values of R\ for positively (left panel) and negatively (right panel) charged 
hadrons as a function of pi . Error bars as in Fig. 2. 



5.1.1 v -dependence 



The first systematic experimental study of the //-dependence of attenuation 
was reported in Ref. [2], where the range 20 < v < 200 GeV was investigated. 
It was shown that the nuclear attenuation decreased with increasing value 
of v, and essentially vanished at v > 50 GeV. The HERMES data are more 
informative, since they are in the region of v where the attenuation becomes 
appreciable, include particle identification, and the statistical precision is much 
better. 



The leftmost columns in Figs. 2, 3, and 5 show that in almost all cases the 
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Fig. 5. Values of R\ for neutral pions as a function of z/, z, Q 2 , and p^. The data 
as a function of z are shown for z > 0.1. Error bars as in Fig. 2. 

attenuation decreases (the value of R\ increases) with increasing values of 
v. (For He this behaviour presumably is present as well, but small compared 
to the uncertainties in the data points.) In the absorption-type models this 
is explained as being due to an increase of the formation length in the rest 
frame of the nucleus at higher v due to Lorentz dilatation, resulting in a larger 
fraction of the hadronization taking place outside the nucleus. In partonic 
models the quark energy loss leads effectively to a shift Az in the argument 
of the fragmentation function, and thus an attenuation that is proportional to 
e/u, with e the quark energy loss. 

For protons R\ increases at higher values of v to well above unity for Kr and 
Xe. Here the following should be realized. The value of (z) is correlated with v, 
e.g., the value of (z) for the lowest z/-bin is about 0.57, whereas for the highest 
u-bin it is 0.35. Since the value of R\ strongly increases with decreasing value 
of z (see the next subsection), a large fraction of the strong increase at high 
v is in fact due to the dependence of R h A on z. Such an effect may play a role 
for other particles, e.g., for K + , too. 



5.1.2 z- dependence 

As can be seen from the second column in Figs. 2, 3, and 5 for all hadron types 
R\ is largely constant with z for He, while it decreases with increasing z for Ne 
and especially for Xe and Kr. In parton energy-loss models this results from 
the strong decrease of the fragmentation function at large z in combination 
with the Az resulting from the energy loss. In absorption-type models the 
overall decrease of R\ with increasing z is assumed to be due to a decrease 
in the formation length in combination with (pre)hadronic absorption. The 
increase of R\ at large z calculated in Ref. [22] is not observed in the data. 
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For the heavier targets R\ rises strongly at low z. Presumably this is due to 
large FSI effects in these nuclei, through which particles of higher energy lose 
energy or get absorbed, generating (other) lower-energy particles. 

As in the case of large values of u, the value of R h A for protons rises above 
unity at small z. This presumably is a result of large rescattering of protons 
and other produced particles with protons in the target. Part of the increase 
is due to the fact that the average value of v decreases considerably with z, 
from about 18 GeV for the lowest bin to about 10 GeV for the highest bin. 
This explains also why even at the smallest z the value of R\ is still lower 
than that for the highest value of v. 

Apart from featuring rather small values of R A , down to almost zero, the z- 
dependence of R\ for antiprotons is special in that it stays constant or even 
decreases slightly at small values of z, where R\ for other particles increases 
strongly for the heavier targets. This probably can be attributed to the fact 
that in final-state interactions the chance that an anti-proton survives, or is 
produced, is relatively small. This would support the idea that the rise of 
R\ for other particles at z < 0.3 is due to FSI effects and suggests that the 
difference in behaviour of anti-protons and the other particles at small values 
of z may be a sensitive check on coupled-channel calculations of FSI effects. 

5.1.3 Q 2 -dependence 

The rightmost column of Figs. 2, 3, and 5 shows for pions a small Q 2 - 
dependence, which is slightly stronger for the heavier nuclei. For kaons and 
(anti-)protons no Q 2 -dependence is discernable. Hence, the attenuation is not 
very sensitive to Q 2 , which means that integrating over Q 2 when studying 
other dependences does not introduce false dependences. 

In the twist-4 energy-loss model of Refs. [20,21] a Q 2 -dependence of R\ of the 
form R ~ a In Q 2 is found, which is consistent with the data. In Ref. [26] the 
calculated Q 2 -dependence is the result of two counteracting processes, which 
results in a rather small Q 2 -dependence that is larger for Kr than for Ne, in 
global agreement with the data. When describing the attenuation purely as 
the result of a modification of the effective fragmentation function [22] , a slight 
increase of R\ with Q 2 is predicted. The deconfinement model [25] predicts 
for all nuclei a slight decrease of R A with Q 2 , which is not supported by the 
data. 

5.1.4 pf -dependence 

Figure 4 and the rightmost column of Fig. 5 show for the heavier nuclei a 
rise of R h A at high p\. Such an effect was first observed by EMC [2] for all 
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charged hadrons taken together, but has been measured now for separate 
identified hadrons. The phenomenon is also known from heavy-ion collisions, 
where it is referred to as the Cronin effect [41]. Compared to experiments 
with heavy ions, the use of a lepton probe has the advantage that initial-state 
interactions do not play a role, except for shadowing effects, which are small 
in the xsj range of the present experiment. The observed rise at high pi is 
attributed to a broadening of the pi distribution. In principle this can result 
from partonic rescattering as well as from hadronic final-state interactions. In 
the next subsection this will be looked at in more detail in the case of pions. 



5.2 Two-dimensional multiplicity ratio for pions 
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Fig. 6. Values of R\ for charged pions for three z ranges. Error bars as in Fig. 2. 
In order to investigate the behaviour of R\ in more detail, the data for 7r H 



20 



and 7r production, which have the best statistical precision and are consistent 
within uncertainties, were combined for a two-dimensional binning, see Figs. 6- 
8. The bins used are 0.2-0.4-0.7-1.2 for z, 6.0-12.0-17.0-23.5 GeV for 
v, and smaller or larger than 0.7 GeV 2 for p 2 . This has the added advantage 
that the correlation between, e.g., the average values of v and z mentioned in 
relation to the one- dimensional distributions, is strongly reduced. 

Figure 6 shows the dependence of R\ on u, Q 2 , and p 2 for three bins in z. The 
left hand column indicates that the dependence on v hardly depends on z. 
The Q 2 -dependence is similar for the various z-bins (Monte-Carlo simulations 
show that the rise of R\ at the lowest Q 2 for the highest z range is due to 
a relatively large contribution of pions coming from p° decay). Therefore, the 
dependence on z is not affected when integrating over Q 2 . The data in the 
rightmost column indicate that the increase of R\ for Kr and Xe at large p 2 is 
smaller for larger z. Such a ^-dependence of the p 2 -dependence was predicted 
in Ref. [26]. The points for the highest z range show a bump at small values 
of pf. This is due to a relatively large contribution of pions coming from p° 
decay in this part of the phase space. 

Figure 7 shows the dependence of R\ on z, Q 2 , and p\ for three bins in v. The 
second and third column indicate that the dependence of R\ on Q 2 and p 2 t 
depends only weakly on which allows one to integrate over these variables 
without introducing spurious correlations. The leftmost column shows for He 
and Ne (and perhaps for Kr) some interesting features in that there seems to 
be a change of the z-dependence with the value of z/, R\ first being about 
constant or even rising slightly with z, and then dropping, the turnover point 
occurring at lower z for the lower z/-bin. This behaviour is studied in more 
detail in the next subsection. 

Figure 8 shows the dependence of R\ on v, Q 2 , and z for two bins in p\. 
The leftmost column shows that at large p 2 t the values of R\ are larger and 
that the dependence on v largely disappears. This is clearly correlated with 
the fact that R h A increases at large p\ (Cronin effect), as discussed in the 
previous subsection, and indicates that this effect is largely independent of v. 
The rightmost column of this figure shows that the Cronin effect disappears 
at high z. This is at least consistent with the idea that the rise of R\ at large 
p 2 (broadening of the ^-distribution) is of partonic origin. In the limit z — > 1 
there is no room for partonic rescattering, because the parton is not allowed 
to have any energy loss (see, e.g., Ref. [26]). In principle, rescattering of the 
produced (pre)hadron could lead to the observed behaviour, too, but since the 
rescattering cross sections are relatively small, the data discussed suggest a 
partonic mechanism. At the same time, this tells that the attenuation in the 
limit z — > 1 is due purely to a hadronic absorption mechanism. 
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5.3 Dependence of R\ on formation length 



Given the ideas of how hadronization proceeds in time, as e.g., in the Lund 
model, the formation length, which depends on both v and z, may be a more 
efficient variable for describing the kinematic dependence of R\ than v and z. 
This idea was recently pursued in Ref . [42] . In order to investigate this, values 
of R\ versus L c for various values of v and z are shown in Fig. 9. Here Eq. 2 
is used with 

f{z) = z^\\-z) (4) 

and k = 1 GeV/fm. This form for f(z) is a convenient parametrization obeying 
the constraints at z — > and z 1, and gives values for L c as a function of 
z closely resembling the ones obtained with the Lund model. 
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A clear correlation can be observed between the values of R\ and L c , with only 
a relatively small residual spread at any fixed value of L c . Evidently most of 
the dependence of R\ on v or z in Figs. 6-7 can be described as a dependence 
on L c , which thus acts as a scaling variable. Upon close inspection it is seen 
that for almost all v values the data points for the lowest two z values (the 
two rightmost in the sequence of same symbols) bend slightly upwards. This 
behaviour is most likely due to large rescattering effects in the lowest z-bin 
(note the relatively strong rise of R\ at low z in Figs. 2 and 3). 

For He the value of R\ rises for the first few points at small L c and then 
becomes constant. For Ne the initial rise extends over more points, and changes 
then to a much more gradual one. For Kr and Xe the change in slope is much 
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Fig. 9. Values of R\{y, z) for He, Kr, Ne, and Xe as a function of the variable L c , 
see Eqs. 2 and 4. The various i/-bins are indicated by different symbols. Within the 
same v bin the z bins used are 0.2 - 0.3 - 0.4 - 0.5 - 0.6 - 0.7 - 0.8 - 0.9 - 1.0, 
and the value of z decreases from left to right. In order not to complicate the figure, 
only statistical error bars are shown. The systematic uncertainty is mainly a scale 
uncertainty of about 3%. 

more gradual, but still noticeable. This is illustrated by the two straight lines 
in the plot for Xe, which represent fits to the data for the ranges L c < 4 fm 
and L c > 5 fm. This suggests the following interpretation: at the larger values 
of L c , which are (much) larger than the size of these nucleJZ] , even if the 
absolute scale of L c may have some uncertainty because the value of k is not 
precisely known, one probably sees a partonic mechanism. The data for Ne 

3 It should be realized that the average distance that a created parton travels 
through a nucleus (assuming it is not absorbed) is only with R the radius of 
that nucleus, because the virtual photon can interact anywhere in a nucleus. Thus, 
even for Kr with a radius of about 5 fm, hadronic mechanisms become small when 
L c > 4 fm. 
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suggest that this mechanism still has a dependence on L c , presumably through 
the underlying values of v and z. The drop of R\ for low values of L c then 
results most probably from hadronic mechanisms. The data on Kr and Xe 
are consistent with this interpretation. At low values of L c there is a strong 
attenuation (stronger in the larger nucleus Xe) due to hadronic absorption (on 
top of the partonic contribution). Since the partonic effect must disappear 
when z — > 1, the value of R\ for L c — > at finite v is purely due to the 
hadronic mechanism. The disappearance of the partonic contribution in this 
limit may be an explanation why the data for the lower values of L c at given 
value of v tend to lie above the data at the same L c with lower i/- value. 
At higher values of L c the influence of absorption becomes smaller, since the 
prehadron is increasingly produced outside the nucleus. Thus above the values 
of L c where the slope changes (below 2 fm for He, and around L c = 2.5, 4, 
and 5 fm for Ne, Kr, and Xe, respectively), one presumably observes mainly 
a partonic mechanism. These tentative conclusions can only be substantiated 
by model calculations that include both partonic and hadronic mechanisms, 
and that give a good description of the measured v and z dependences of R\. 
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Fig. 11. The dependence of the parameters a and [3 on the value of L c for 
0.3 < z < 0.7. 



5.4 A- dependence 



The /L-dependence of R\ was fitted in a phenomenological way by using as 
the fit function 

R h A (A)=exp[-P(A/100n (5) 

The scale factor of 100 is introduced to reduce the correlation between a and 
(3. At the same time the value of f3 is representative now of the attenuation of 
a nucleus with A = 100. For moderate values of (3A a this formula is equivalent 
to 1 — R\ = /3A a , a parametrization that has been used before [1], but the 
present form has the advantage that R\ does not become unphysical for very 
large values of A. 

However, since R\ is the ratio of the multiplicity in nucleus A vs. the mul- 
tiplicity in the deuteron, the fit function of Eq. 5 is internally inconsistent, 
unless one assumes no attenuation in the deuteron. Assuming that the at- 
tenuation depends on the (average) density times the radius of the nucleus, 
in Ref. [43] it is found that the attenuation in a deuteron can be described 
with an effective A-value of about 0.6. The influence of this on the value of 
a obtained using Eq. 5 is fairly small (about -0.07). Therefore, all results to 
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Fig. 12. Values of the parameter a for pions, kaons, and (anti-)protons in three v 
ranges (left six panels), and in three z ranges (right six panels). 

be discussed in the following are based on Eq. 5, but in drawing conclusions 
it should be realized that the real value of a is probably slightly lower. More 
information on the performed fits with different type of fitting functions can 
be found in [44]. 

A second remark concerns what A-dependence is expected from different types 
of models for the attenuation. As already mentioned in section 2, the attenu- 
ation in the parton energy-loss model of Refs. [20,21] is given as 1 — Ra ~ L 2 ~ 
A 2 / 3 , whereas in (Glauber type) absorption models it is often presumed that 
1 — Ra ~ L ~ A 1 / 3 . However, these estimates are too simple. Taking realistic 
matter distributions of nuclei yields effective values of L 2 that are proportional 
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to A 0,74 [43]. Furthermore it was demonstrated already in Ref. [25] that the 
inclusion of a distribution for the formation length in absorption calculations 
increases the exponent above the value of 1/3, yielding values for 1 — that 
approximately follow an A 2 / 3 pattern for large values of L c , if the nucleus is 
described as a sphere with constant density. Using realistic matter distribu- 
tions the authors of Ref. [43] find the exponent of A to be 0.40, 0.54, and 0.60 
when L c = 0, 2, and 4 fm, respectively. 

In the fits only the statistical uncertainties in the values of R\ have been used, 
as the systematic uncertainty in the values of R\ is largely a scale uncertainty. 
The influence of the latter on the value of a was found to be about 0.05. 

5.4-1 Pions 

First the A-dependence of pion production was studied as a function of v and 
z, using the data binned in both z and v. The results for a as a function of z 
for the various i^-bins are shown in Fig. 10. The values of (3 (not shown) reflect 
the global behaviour of R\ (more attenuation at higher z and lower v) that 
is visible from earlier figures. There is an increase of a from about 0.5 to 0.6 
with z and possibly at higher z a slight increase with v. However, in the latter 
region the results may be influenced by contributions from the decay of the p°. 
The results for the highest z-bin do not follow the pattern of the lower z-bins. 
This is due to the behaviour of R h A for especially He and Ne, seen in Fig. 7. 
As mentioned before, this could result from the formation length becoming so 
small that even in He and Ne hadronic mechanisms become important. 

Given these results for the separate v- and ^-dependences it has also been 
investigated how a and (3 depend on the value of L c . For that purpose the 
data were binned in five L c -bins, using the values of v and z of each event in 
Eqs. 2 and 4. In order to avoid possible contributions of the p° at high z and 
large rescattering effects at low z, only data with 0.3 < z < 0.7 were used. The 
resultant values of a and (3 are shown in Fig. 11. The behaviour of (3 reflects 
the smaller attenuation at higher L c already visible in Fig. 9. Consistent with 
what was found above, the value of a decreases from about 0.6 at small L c to 
less than 0.5 at large L c . 

If, as was argued in the previous subsection, partonic effects are most promi- 
nent at large values of L c , this would indicate that the A-dependence of that 
mechanism has a value of a well below the value of 2/3 given in Refs. [20,21] 
and [22]. 

The decrease of a with L c is not easily explained. As demonstrated in Refs. [25,43] 
a pure absorption mechanism would yield a value of a that increases with L c . 
Possibly such an increase is more than compensated for by an increasing influ- 
ence of the partonic mechanism. However, then that mechamism should have 
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an A-dependence with a rather small value of a. But, when two (or possibly 
even more) mechanisms contribute, the A-dependence probably becomes more 
complicated than can be described by a single exponential as in Eq. 5. For 
that reason comparisons between data and model calculations should be done 
on the level of the multiplcity ratio R\. 

5.4.2 Other particles 

For the other particles the statistical precision is too low for a meaningful 
two-dimensional binning and study of the A-dependence, so there the A- 
dependence is presented for i/-bins of 6.0 — 12.0 — 17.0 — 23.5 GeV (integrated 
over all z), and for z-bins of 0.2 — 0.4 — 0.7 — 1.2 (integrated over all v). For 
comparison the values for ir + and ir~ in the same bins are given as well. The 
results are shown in Fig. 12. The only significant feature is the different be- 
haviour seen for K + particles at high v. As mentioned before, this may be due 
to large rescattering effects. The behaviour for K~ is different from the one 
for K + , and more like that of 7r _ , but with smaller values for a. 



6 Summary and conclusions 

Data for the multiplicity ratio R\ of hadron production in semi-inclusive deep- 
inelastic scattering of 27.6 GeV electrons and positrons from helium, neon, 
krypton, and xenon nuclear targets relative to deuterium were obtained for 
identified n + , n~ , ir°, K + , K~ , p, and p particles as a function of the virtual- 
photon energy u, the fraction z of the energy transferred to the hadron, the 
photon virtuality Q 2 , and the hadron transverse momentum squared p 2 . For all 
particles the dependence of R\ on these variables is presented and discussed. 

The most prominent features of the data are an increased attenuation (de- 
crease of R\ below unity) with increasing value of the mass number A of the 
nucleus and the attenuation becoming smaller (larger) with increasing values 
of v (z), R\ dropping below 0.5 for xenon in some kinematic regions. At low 
values of z, especially for heavier targets and for protons and K + particles, a 
strong rise of R\, even to above unity, is observed. Presumably this is due to 
hadronic rescattering, where a higher energy particle through nuclear reactions 
produces one or more lower-energy particles. 

The value of R\ increases slightly with Q 2 , at least for pions, and is almost 
independent of pf, except at large values of p 2 , where R h A increases strongly. 
The latter is thought to result from p\ broadening due to partonic rescattering 
(Cronin effect). This effect was seen to disappear for z 1, in accordance with 
the picture that in that limit no rescattering is possible, since rescattering of 
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the struck parton implies an energy loss. 

By combining the data for 7r + and ir~ , the dependence of R\ on two of the 
variables v, z, Q 2 , and p 2 together was investigated. The dependence on Q 2 
depends weakly but noticeably on the value of u, but practically not on that 
of z. The dependence on p 2 hardly depends on v and z, except for the dis- 
apperance of the rise at large p 2 at z — > 1 mentioned above. However, the 
dependences on v and z are related. It was found that most of the depen- 
dence on v and z can be incorporated in a dependence on the combination 
L c = z 0,35 (l — z)v/k, where k is the string tension in string models, which 
thus acts as a scaling variable. Since this function is close to the one given in 
the Lund model for the average formation length of a particle, by inspecting 
the value of R\ vs. L c for the four nuclei, regions can tentatively be identified, 
where hadronic (absorption) plus partonic mechanisms are important, and a 
region at higher L c where only or mainly partonic mechanisms play a role. 

A fit of the A-dependence of the values of R\ for pions measured for the various 
nuclei of the form R\ = exp[— /3(A/100) Q ] yields values of a from about 0.6 
to 0.5, depending on the value of L c . Similar values are found for the other 
particles. These values of a are well below the values resulting from models in 
which the attenuation depends on the square of the distance a parton travels 
through the nucleus. 

In total a very extensive data set to guide modeling hadronization in nu- 
clear matter has been collected. A full theoretical description of hadronization 
in nuclei in one consistent framework, including partonic and hadronic (ab- 
sorption plus rescattering) mechanisms is badly needed. Clearly it will be a 
challenge for any theoretical model that is developed to describe these data 
for the various hadrons and nuclei as a function of all kinematic variables, but 
if successful, this combination of data and theoretical interpretation will con- 
tribute essentially to the understanding of non-perturbative QCD at normal, 
and thence higher densities. 
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A Appendix 

In the next table the average values of the kinematic variables that were inte- 
grated over when showing the dependences on is, z, Q 2 and pf in Figs. 2, 3, 4 
are given for the case of pion production on krypton. These values hardly 
depend on the target used. 
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Table A.l 

Average values of v, Q 2 , z, and p 2 for pions produced on krypton. 



v range 


(u) (GeV) 


(*> 


(Q 2 ) (GeV 2 ) 


(p 2 ) (GeV 2 ) 


4.0 - 6.0 


5.269 


0.583 


1.889 


0.245 


6.0 - 8.0 


7.151 


0.501 


2.029 


0.164 


8.0 - 10.0 


9.092 


0.448 


2.203 


0.141 


10.0 - 12.0 


11.030 


0.421 


2.352 


0.147 


12.0 - 14.0 


13.009 


0.406 


2.506 


0.171 


14.0 - 16.0 


14.991 


0.392 


2.595 


0.199 


16.0 - 18.0 


16.980 


0.370 


2.584 


0.234 


18.0 - 20.0 


18.964 


0.349 


2.475 


0.275 


20.0 - 23.5 


21.540 


0.326 


2.164 


0.334 


z range 


<2> 


(u) (GeV) 


(Q 2 ) (GeV 2 ) 


(p 2 ) (GeV 2 ) 


0.1 - 0.2 


0.155 


18.358 


2.463 


0.104 


0.2 - 0.3 


0.247 


15.953 


2.623 


0.164 


0.3 - 0.4 


0.345 


14.796 


2.655 


0.219 


0.4 - 0.5 


0.446 


14.315 


2.635 


0.260 


0.5 - 0.6 


0.546 


13.963 


2.585 


0.284 


0.6 - 0.7 


0.646 


13.637 


2.500 


0.286 


0.7 - 0.8 


0.746 


12.824 


2.392 


0.246 


0.8 - 0.9 


0.845 


11.841 


2.252 


0.192 


0.9 - 1.3 


0.952 


10.742 


2.258 


0.156 


Q 2 range 


(Q 2 ) {GeV 2 ) 


(v) (GeV) 


<*) 


(p 2 ) (GeV 2 ) 


1.0 - 1.5 


1.360 


14.790 


0.399 


0.217 


1.5 - 2.0 


1.911 


14.646 


0.395 


0.212 


2.0 - 3.0 


2.667 


14.874 


0.389 


0.213 


3.0 - 4.0 


3.746 


15.062 


0.382 


0.213 


4.0 - 5.0 


4.775 


15.100 


0.380 


0.211 


5.0 - 6.0 


5.798 


15.065 


0.380 


0.208 


6.0 - 8.0 


7.161 


14.819 


0.380 


0.203 


8.0 - 25.0 


9.735 


14.583 


0.386 


0.192 


p 2 range 


(p 2 ) {GeV 2 ) 


(u) (GeV) 




(Q 2 ) (GeV 2 ) 


0.00 - 0.05 


0.023 


13.059 


0.379 


2.550 


0.05 - 0.10 


0.073 


14.107 


0.390 


2.611 


0.10 - 0.30 


0.182 


15.276 


0.385 


2.599 


0.30 - 0.50 


0.384 


16.172 


0.392 


2.616 


0.50 - 0.70 


0.585 


16.669 


0.424 


2.635 


0.70 - 0.90 


0.787 


16.912 


0.456 


2.627 


0.90 - 1.10 


0.987 


17.211 


0.483 


2.620 


1.10 - 1.40 


1.225 


17.471 


0.505 


2.648 


1.40 - 1.95 


1.606 


17.788 


0.535 


2.615 


1.95 - 5.00 


2.418 


18.487 


0.571 


2.625 
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